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INTRODUCTION 

Finite  state  space  semi-Markov  processes  find  application  in  many  areas 


Including  queueing,  reliability,  and  clinical  trials;  (cf  Cinlar  [1975],  Weiss 


and  Zelen  [1965],  Cox  [1985]).  Often  Interest  in  the  application  centers  on 
the  distribution  of  a  fl  rst -passage  time  to  a  state  or  a  set  of  states 
representing  for  example  the  lifetime  of  a  system  or  the  end  of  a  busy  period 
of  a  server. 

In  this  paper  we  consider  the  problem  of  the  estimation  of  the  log 
survivor  function  of  a  first  passage  time  for  a  fixed  time  t.  We  suppose  that 
observational  data  are  known  about  the  semi  Markov  process  in  question. 
Although  the  approach  and  results  are  given  concretely  for  one  particular 
semi -Markov  process,  they  apply  more  widely. 


To  be  specific,  consider  N  individuals.  Let  X^(i)  be  the  state  of  1th 


individual  at  time  t.  We  will  assume  {X^(i);t>0}  are  ind(  pendent  identically 


distributed  processes  having  the  same  probability  law  as  {X^;t>0) 


The 


process  {X^:t>0}  is  a  semi-Markov  process  with  three  states  {0,1, 2).  State  0 

The  sojourn  time  in  state  1  has  a  distribution 
Upon  leaving  state  1,  the  process  transitions  to  state  0  with 


is  absorbing  and  Xq=1. 


function  F 


i 


probability  9  and  to  state  2  with  probability  1-9 
transitions  to  state  1  with  probability  1. 


From  state  2  the  process 
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The  problem  Is  to  estimate  the  logarithm  of  P{D>t),  for  fixed  t.  from  data 
obtained  by  observing  the  N  individuals.  Three  possible  estimation  procedures 
will  be  considered.  The  three  estimators  use  different  amounts  of  information 
concerning  the  process.  One  procedure  uses  only  the  observed  entrance  time  to 
state  0  for  the  N  individuals  to  estimate  the  empirical  distribution  function 
of  D.  Another  procedure  makes  parametric  assumptions  concerning  the  sojourn 
time  distributions  Fj  and  uses  maximum  likelihood  to  estimate  the  probability. 
A  third  approach  uses  an  exponential  approximation  to  P(D>t}  and  empirical 
distributions  to  estimate  F^. 

In  Section  2.  the  estimation  procedure  based  on  an  exponential 
approximation  »o  P{D>t)  is  described  and  asymptotic  results  are  obtained  for 
it. 

In  Section  3  the  other  two  estimation  procedures  are  described  and  the 
results  of  a  simulation  experiment  are  given.  The  experiment  was  conducted  to 
compare  the  performance  of  confidence  intervals  for  the  three  estimators  for 
small  and  moderate  numbers  of  individuals.  It  is  found  that  the  procedure 
based  on  the  maximum  likelihood  is  the  most  efficient  if  the  underlying  model 
is  correct.  However,  maximum  likelihood  confidence  Interval  coverage  is 
sensitive  to  incorrect  model  assumptions.  The  nonparametrlc  confidence 
intervals  based  on  the  exponential  approximation  have  the  correct  coverage  as 
long  as  the  time  t  is  not  too  small.  The  conservative  binomial  confidence 
intervals  based  only  on  the  observed  times  of  absorption  are  the  least 
efficient  and  tend  to  overcover. 

2.  AN  ASYMPTOTIC  RENEWAL  ESTIMATE 

The  probability  P{D>t)  for  the  model  described  in  the  Introduction 


2 


satisfies  the  equation 


P{D>t)  =  g(t)  +  (l-«)  /q  (Fj^FgXds)  P{D>t  s} 


(2.1) 


g(t)  -  Fj(t)  4  (i-«)  Fj((ls)  Fgit-s). 


(2.2) 


This  is  a  renewal  equation  with  defective  inter -renewal  di.stribution 


L(t)  (l-«)(Fj*F2)(t) 


(2.3) 


Following  Feller  [1966]  page  326,  assuae  there  exists  x  which  satisfies 
the  following  equation 


/  e""  I,(dt)  =  (1-S)  ♦  (X)  ♦  (K)  =  1 

0 


(2.4) 


where 


♦  ({)  =  /  e'"  F  (dt) 

0 


(2.5) 


If  arithaetic,  then  under  certain  integrabl 1 Ity  conditions 

lim  e"*^  P{D>t)  =  ^  (2.6) 


3 


U*JI  Wjnj)  15 


where 


and 


W 


Oft  if  o 

/i  =  /  se  L(ds) 
0 


(2.7) 


b  =  /*  e'*''  g(3)ds  =  1  ♦,(«).  (2.8) 

0 

Following  an  approach  that  was  found  useful  in  another  context,  (Gaver 
and  Jacobs  ([1986]),  an  estlaator  for  P{D>t)  can  be  obtained  by  using 

estinates  of  x,  b  and  fJ  In  (2.6).  Specifically,  let  (S^(l),n=l .  M^}  be 

the  collection  of  sojourn  tiwes  in  state  i  for  the  N  individuals  and  let  R  be 
the  nuaber  of  transitions  fruM  state  1  to  state  2  for  all  the  N  individuals. 


p„, , .  _  .„d 


;  ,t,  >  ?  e'V" 

*  k=l 


1  -  i,a 


An  estlaator  of  x  is  the  solution,  x,  to  the  equation 


(2.9) 


(1-®)  ♦j(e)®2(t)  =  1  (2.10) 

Notice  that  the  left  hand  side  of  equation  (3.10)  is  nondecreasing  in  C  and  at 
(  ^  0  equals  (1  ®)  <  1.  Hence  (2.10)  has  a  unique  solution  which  can  be  found 
nuaerical ly . 
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The  estliiators  for  the  other  quantities  In  (2.6)  are 


«i  ”2 


^  ^  exp{/c(S.  (1)+S,(2)))  (S.  (1)  +  S.(2))  (2.11) 


«l”2 


k '  *  '  j ' 


J' 


k=l  J-l 


0 


b  -  TT  ♦  j  )  . 

K 


(2.12) 


The  asymptotic  exponential  estimator  for  P{D>t)  is 


P^{D>t}  =  ^e'^V 
A* 


(2.13) 


2 . 1  Large  Sample  Properties  of  k  and  fn  P^(J?>U 


Let 


f(e.a)  =  1  -  (l-a)#j(C)02({)  2  1  (l-a)h(C). 


Expand  h(x)  In  Taylor  series  about  the  solution  k  of  (2.4).  Since  f(x.0)  =  0, 


0  =  f(x.0)  =  1  (1-0)  [h(K)  t  (ic-x)hM«)  +  h”(ax)]  (2.14) 


=  f(K,0)  +  (0-0)h(K)  (1-0)(K-/C)[h'  («)  +  i(K-K)h'  ■  (cUC)l  . 


Thus 


f(K,0)  ♦  (0-0)  h(K) 


K  -K  = 


(2.15) 


where 


B  -  (l-«)th*('‘)  +  h“(ax)]. 


(2.16) 


and  a  is  a  constant  and  «  is  the  probability  of  a  transition  from  state  1  to 
state  0.  By  the  strong  law  of  large  nuwbers.  as  f(#c,e)-»0  and  9-0-*0. 
Further,  as  Nhp»,  with  probability  one 


and 


h'(«)  -*  (l-«)  se  ‘F  )(d8) 

0  ^ 


h(/{)  -  ♦j(K)  ♦gl*)- 


(2.17) 


(2.18) 


Therefore,  k-«-«0  as  NM»,  with  probability  1  so  the  estimator  k  is  consistent. 

To  obtain  a  central  limit  theorem  for  «.  the  following  moment  condition 
is  needed.  Assume 


r  F  *F  tds)  <  «. 

0 


(2.19) 


Since  f(ic,8)  =  0, 


VTT  [  f(«,«)]  =  ^  [  f(w.®)  f(«.«)) 


N  (l-«)  (♦j(k)^2(«)  ♦i(«K2<'‘)J 


yH  (1-6)  ((♦j(«)  -  ♦j(K)]  ♦gdt)  ♦  ♦j(«)[62('')^2^'‘*^^- 

6 


,S' 

.V 


Thus,  (2.19)  and  the  central  limit  theorem  imply  that  y~N  (f('t.9)l  is 
asymptotically  normal  as  N-io*  with  mean  0  and  variance 

K  T  K  T 

Var^  =  (]-«)  9{[-i  Varfe  Var(e  ^]}  (2.21) 

where  T.  is  a  random  variable  with  distribution  F.. 

1  1 

Further  the  central  limit  theorem  implies  that  as  N-*».  (®-®)  is 

2  2  3 

asymptotically  normal  with  mean  0  and  variance  =  9  (1-9)  *  9  (1-9). 

Hence,  as  N-»»,  {k  k)  is  asymptotically  normal  with  mean  0  and  variance 

-  (V^  +  (♦^(K)  Vj^)[(1-9)  E[(Tj  +  T2)e'‘^'’'l''‘'2^]]'^  (2.22) 

A  similar  argument  shows  that  P^{U>t}  is  also  asymptotically  normally 
distributed  as 

3.  SIMULATION  RESULTS 

A  simulation  experiment  was  carried  out  to  investigate  small  and  moderate 
sample  size  behavior  of  confidence  Interval  procedures  for  three  estimators 
of  fn  P{D>t).  The  following  three  estimators  use  different  amounts  of 

Information  about  the  semi-Markov  processes  (X^(l):  t>0),  i=l .  N.  for  the 

N  individuals 

3.1  The  Binomial  Estimator 

An  estimator  which  uses  only  the  observed  entrance  times  into  state  0  for 

the  N  individuals,  d,,d  , .,.,d„,  is 

1  c  N 
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(3.1) 


where 


1  if  d^>t. 


0  otherwise. 


The  distribution  of  the  estiaator  P  {D>t}  is  binonial  and  a  binomial 

O 


confidence  interval  [L  .1)  ]  can  be  constructed  for  P(D>t)  cf.  Larson  [1982]  p 

B  B 


397.  The  binomial  confidence  Interval  for  fn  P{D>t}  is  [<n  L  ,fn  U  ]. 

B  B 


3.2  The  Asymptotic  Exponential  Estimator. 


The  asymptotic  exponential  estimator  for  P{D>t}  is  given  by  P^{D>t)  of 


(2,13).  The  estimator  for  the  fn  P{D>t)  is  fn  P.{D>t)  and  is  asymptotically 

A 


normal  as  However,  since  the  simulation  experiment  is  for  small  to 


moderate  numbers  of  individuals,  two  nonparametr ic  methods  are  used  to  obtain 


confidence  Intervals, 


The  jackkfiife  is  a  procedure  originally  introduced  by  Quenoullle  [1956] 


for  bias  reduction,  and  adapted  by  Tukey  [1958]  to  obtain  approximate 


confidence  intervals.  Suppose  Interest  is  on  a  parameter  p  that  is  estimated 


by  p  u.sing  a  complex  calculation  from  data.  The  idea  is  that  of  assessing 


variability  by  recomputing/!  after  removing  independent  subgroups  of  data  and 


then  using  the  recomputed  p  values  to  estimate  a  variance  which  is  In  turn 


applied  to  state  a  two-sided  confidence  Interval  that  contains  the  true  p  with 


specified  confidence:  see  Efron  [1982J  and  his  more  recent  work,  or  Hosteller 
and  Tukey  [1977]  for  more  details.  For  simulations  Involving  N  20 
individuals,  the  jth  subgroup  consists  of  all  data  corresponding  to  the  jth 
individual.  For  simulations  involving  N=50  Individuals  the  first  subgroup 
consists  of  all  data  corresponding  to  the  first  five  individuals  the  second 

subgroup  of  all  data  corresponding  to  the  second  five  itidi viduals .  etc.  Some 

simulations  for  N=50  individuals  were  run  with  the  jth  subgroup  consisting  of 
the  data  for  the  jth  Individual  alone;  according  to  some  work  these  should  be 
an  improvement  on  those  with  5  in  a  subgroup.  However,  the  resulting 
confidence  Intervals  differed  little  from  those  obtained  by  leaving  out  5 

individuals  at  a  time. 

The  bootstrap  is  an  alternate  method  that  may  be  used  to  assess  uncertainty 
of  estimators:  of.  Efron  [1982]  and  his  more  recent  work.  In  the  simulations 
reported  here  a  bootstrap  replication  is  generated  as  follows.  Let  be  the 
number  of  transitions  from  state  1  to  state  2  for  individual  n.  Let 
be  the  collection  of  all  sojourn  times  in  state  i  for  all  individuals.  To 
generate  bootstrap  data  for  one  individual,  one  observation  is  drawn  at  random 
with  replacement  from  call  the  observation  r;  r«^l  observations  are  drawn 

at  random  with  replacement  from  the  collection  of  state-1  sojourn  times 
{S,  (1));  r  observations  are  drawn  at  random  with  replacement  from  the 

K 

collection  of  state  2  sojourn  times  {Sj^(2)).  Bootstrap  data  is  generated  for 
N  individuals  and  the  estimate  fn  P.{D>t)  is  computed  using  this  data.  This 
completes  one  bootstrap  leplication.  B  =  100  bootstrap  replications  are  done. 
The  B  =  100  bootstrap  estimates  of  fn  P^{D>t)  are  ordered.  A  100  x  (1-a)  % 
confidence  interval  is  constructed  using  the  a/2  and  1 -a/2  quantiles  of  the 
bootstrap  estimates 
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In  what  follows,  we  assume  the  sojourn  time  in  state  i  has  an  exponential 


distribution  with 
assumptions  is 


1 


The  log-likelihood  function  under  these 


L 


N 

=  Z 
n=l 


[R^<n  (1-®)  <n  e  +  R^fn  + 


(1+R^)«n 


(3.2) 


where  R  is  the  number  of  transitions  from  state  1  to  state  2  for  individual  n 
n 

and  T^(l)  is  the  total  time  spent  in  state  1  by  individual  n.  The  maximum 
likelihood  estimators  are 


N 

Where  R  =  Z  R  ,  T(l) 
n 

n=l 


N^R 

*  N+R 

'  Tfl) 

^2  "  fT2) 

N 

Z  T  (1 ) ,  and  T(2) 
n 

n=l 


N 

Z  T  (2)  . 
n 

n^l 


(3.3) 

(3.4) 

(3.5) 


Let  D  be  the  entrance  time  to  state  0  for  an  individual;  then 


Ag+Pg  A  4p  9p 

S(t)  s  P{D>t)  =  ^  exp{A  t)  -  exp{A  t}}  r— ^  (3.6) 

2  2  Aj  1  Aj  A^ 
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where  ^j.*2  roots  of  tht;  equation 

•PjPg  "  y(Pi*p2’  +  =  0.  (.1.7) 

The  waxiaun  likelihood  estiaator  for  P{D>t}  can  be  obtained  from  (3.6)  and 
(3.7)  by  replacing  .  p^  and  0  by  their  aaxiaua  likelihood  estiaators  and 
taking  the  logaritha  of  the  resulting  estiaator  of  F'{n>t)  Tfie  distribution 
of  the  resulting  estiaator  is  asyaptotical ly  noraal  as  N-^.  Since  the  aaxiaua 
likelihood  estiaators  ^nd  p^  are  uncorrelated,  the  asyaptotlc  variance 

is 


Var[#n  P(D>t}] 


S(t)*  ' 


Var(p2 


2 

)  J.  (3.8) 


Confidence  intervals  are  constructed  using  the  asyaptotic  noraal ity  of  the 
aaxiaua  likelihood  estiaator  and  the  asyaptotic  variance  with  paraaeter  values 
being  replaced  by  their  corresponding  estiaates.  Confidence  intervals  based 
on  the  aaxiaua  likelihood  ratio  test  aay  also  be  constructed  [cf.  Cox  and 
Hinkley  (1974)  pp.  343].  While  these  intervals  are  asyaptot ii:ally  equivalent 
to  the  noraal  confidence  intervals  that  are  computed,  they  aay  be  better  for 
saaple  sizes. 


3.4  ihtf  slBulaUQii  Expgriient 

All  siaulations  are  carried  out  on  an  IBM  3033AP  coaputer  at  the  Naval 


Postgraduate  School  using  the  LLRANDONII  randoa  number  generating  package;  [cf 


Lewis  and  Uribe  (1981)].  Soae  details  of  the  simulation  are  given  below.  A 
more  complete  account  can  be  found  in  Kin  [1987].  In  each  replication 
estimates  and  confidence  intervals  are  computed  for  the  logarithm  of  P([)>t) 
using  the  procedures  described  in  subsections  3. 1-3.3.  Each  simulation  is 
replicated  300  times.  For  each  procedure,  the  number  of  confidence  intervals 
covering  the  true  value  of  fn  P{D>t}  is  recorded.  Also  recorded  are  the 
number  of  intervals  that  are  too  low  (true  tn  P{D>t}  >  U,  the  upper  endpoint 
of  the  interval)  and  too  high  (true  fn  P{D>t}  <  L,  the  lower  endpoint  of  the 
interval).  The  average  length  of  the  confidence  Intervals  is  also  computed. 

Tables  1-8  present  results  from  simulations  for  N=20  and  N=50 
individuals.  The  first  row  in  each  table  lists  the  confidence  interval 
procedures:  the  binomial  confidence  interval  for  logarithm  of  estimate  (3.1) 
(BIN):  the  maximum  likelihood  Interval  (NLE);  the  jackknife  confidence 

interval  for  the  asymptotic  estimate,  (JK);  and  the  bootstrap  confidence 
interval  for  the  asymptotic  estimate,  (BT).  Nominal  80%  and  90%  confidence 
interval  for  fn  P(D>t)  are  computed  using  each  procedure  for  each  time 
*■  =  0.5,  1.0,  1.5,  2,  3,  4.  Confidence  limits  that  are  greater  than  0  are  set 
equal  to  0. 

Tables  1-4  present  results  for  the  case  in  which  the  sojourn  times  in 
state  1  have  an  exponential  distribution  with  =  1  and  the  sojourn  times  in 
state  2  have  an  exponential  distribution  with  mean  =  0.1.  The  values  of 

P{D>t)  for  the  model  for  t  =  0.5,  1.0,  1.5,  2,  3,  4  are  given  in  Appendix  A. 

Tables  1-2  show  results  for  simulations  with  N=20  individuals. 

Table  1  shows  coverage  results.  Under  each  procedure  is  given  the  number  of 
intervals  that  cover  the  true  value  (C),  the  number  of  intervals  that  are  too 
low  (L)  and  the  number  of  Intervals  that  are  too  high  (H).  In  parentheses 


next  to  these  numbers  are  given  the  corresponding  fractions. 


If  the  true 


confidence  level  of  a  confidence  Interval  procedure  is  80!fe,  (respectively 
90!li)  .  then  a  95%  confidence  interval  for  the  fraction  of  the  300  replications 

that  cover  the  true  value  is  0.80  ±  1.96  (.2)(.8)  =  [.755,  .845], 

(respectively  .9  ±  1.96  ^  (.!)(  9)  =  [.866,  .934]). 

Table  2  shows  the  average  length  of  confidence  Intervals  for  fn  P(D>t} 
for  each  t  and  each  procedure.  In  parentheses  under  the  average  length  is  the 
estimated  standard  error  of  the  average  length. 

The  conservative  binomial  procedure  tends  to  overcover  and  has  the 
largest  average  confidence  interval  length.  The  maximum  likelihood  procedure 
has  the  correct  coverage  and  smallest  average  length.  The  confidence 
intervals  based  on  the  asymptotic  renewal  estimate  show  the  correct  coverage 
for  the  times  larger  than  t  =  0.5  and  1.0.  When  the  Jackknife  and  bootstrap 
Intervals  do  not  cover  they  tend  to  be  too  high.  The  average  length  of  the 
Jackknife  and  bootstrap  intervals  is  between  those  of  the  binomial  and  maximum 
likelihood  intervals.  The  average  length  of  the  Jackknife  intervals  is 
greater  than  those  for  the  bootstrap,  suggesting  that  the  latter  is  preferred: 
however  computational  cost  may  be  higher. 

Tables  3  a/id  4  present  coverage  results  and  average  confidence  interval 
lengths  for  a  simulation  experiment  with  N=50  individuals.  The  coverage 
results  are  similar  to  those  of  Table  1.  As  anticipated,  the  average 
confidence  interval  lengths  tend  to  be  smaller  for  N  =  50  individuals  than  for 
N=20  individuals.  The  average  bootstrap  interval  length  is  closer  to  that  of 
the  MLE  interval  for  times  larger  than  1  for  N=50  individuals  than  for  N=20 
individuals . 

Tables  5-8  report  results  of  a  simulation  experiment  in  which  the  sojourn 

13 


tiKes  in  states  1  and  2  have  the  following  distributions.  The  sojourn  tine  In 


state  2  has  an  exponential  distribution  with  mean  —  =  0.1  as  before.  The 

^2 

sojourn  tiae  In  state  1  Is  the  sua  of  two  independent  exponential  randoa 
variables  each  having  aean  j.  Thus,  the  aean  sojourn  tiae  In  each  state  Is 
the  saae  as  for  the  aodel  used  In  Tables  1-4.  Values  of  P{D>t}  for  various 
values  of  t  for  this  aodel  are  given  in  Appendix  B.  Estiaators  and  confidence 
Intervals  for  tn  P{D>  )  are  coaputed  as  before.  In  particular  the  aaximua 
likelihood  estiaators  assuae  that  the  sojourn  tiaes  in  state  1  have  an 
exponential  distribution  with  aean  1.  The  nuaber  of  individuals  is  N-20  for 
the  siaulation  in  Tables  5-6.  There  are  N-50  individuals  in  the  simulations 
of  Tables  7-8.  The  coverage  results  contained  in  Tables  5  and  7  show  that  the 
aaximua  likelihood  confidence  Intervals  are  sensitive  to  underlying  model 
assumptions.  Incorrect  model  assumptions  lead  to  both  overcoverage  and 
undercoverage.  The  confidence  intervals  based  on  the  exponent  ial 

approximation  have  the  correct  coverage  for  the  times  larger  than  0.5.  The 
average  length  of  the  Jackknife  and  bootstrap  intervals  is  smaller  than  those 
for  the  binomial.  The  average  length  of  the  bootstrap  Interval  is  smaller 


than  those  for  the  Jackknife. 


4.  CONCLUSIONS 

The  simulation  results  indicate  that  the  maximum  likelihood  estimator  is 
most  efficient  when  the  correct  model  is  used.  As  has  beijn  noted  elsewhere, 
the  MLE  estimator  is  rather  sensitive  to  Incorrect  aodel  assumptions.  The 
simulation  results  for  the  asymptotic  exponential  estimator  are  somewhat 
surprising  in  that  it  appears  that  the  time  t  does  not  have  to  be  very  large 
in  order  for  the  asymptotic  exponential  estimator  to  work  well.  Both  the 


nonparaaetric  bootstrap  and  jackknife  confidence  Intervals  for  the  asymptotic 
exponential  estimator  have  the  correct  coverage  for  moderately  large  t.  As 
implemented  the  bootstrap  Intervals  tend  to  have  a  smaller  average  length  than 
the  jackknife.  However,  neither  a  complete  Jackknife  nor  a  bootstrap  with 
more  than  100  replications  has  been  done  due  to  the  massive  amount  of 
computation  involved.  As  expected,  the  binomial  estimator  produces 
conservative  confidence  Intervals  which  tend  to  overcover.  However,  the 
binomial  confidence  Intervals  are  appropriate  for  all  times  and  all 
semi -Markov  models. 
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Appendix  A 


Values  of  P{D>t}  for  the  seml-Karkov  Model 
with  exponential  sojourn  time  in  state  1 

e=0.5 


Ppl. 


P2=10. 


Appendix  B 


Values  of  P{D>t}  for  the  semi-Markov  Model 
with  nonexponential  sojourn  time  in  state  1 
p^=2,  P2=10.  P3=2.  0=0.5 


time  t 

PID>tl 

0.5 

.8652 

1.0 

.6743 

1.5 

.5157 

2.0 

.3930 

3.0 

.2278 

4.0 

.1320 

Table  1 


Coverage  Results  for  Confidence  Intervals  for  #n  P{D>t} 
Exponential  Sojourn  Times 
N=20.  0=0.5.  pj=1.0.  p2=10.0 


Time 

Conf idence 
Level 

Coverage 

BIN 

MLE 

JK 

BT 

0.5 

SOX 

H 

14{.05) 

24(.08) 

106{.35) 

140(.47) 

C 

272{.91) 

251(.S4) 

174(.5S) 

160(.53) 

L 

14^05) 

25(.08) 

20^07) 

0(0) 

90X 

H 

4(01) 

14(.05) 

7S{ .26) 

99(.33) 

C 

2S2(.94) 

276(.92) 

211(.70) 

201^67) 

L 

14^05) 

10(.03) 

11(.04) 

1.0 

SOX 

H 

17(.06) 

24(.0S) 

72(.24) 

66(.22) 

C 

250(,S3) 

250(.S3) 

210(.70) 

224(.75) 

L 

33^11) 

26(.09) 

1S(.06) 

10{.03) 

90X 

H 

7(,02) 

14(.05) 

46( . 15) 

33{.ll) 

C 

275(.92) 

277 (.92) 

245(.S2) 

254(.86) 

L 

18(.06) 

9(.03) 

9(.03) 

8(.03 

1.5 

SOX 

H 

15(.05) 

24(.0S) 

4S(.16) 

40{. 13) 

C 

265(.SS) 

250(.S3) 

247(.82) 

236( .79) 

L 

20(.07) 

26(.09) 

5{.02) 

24(.08) 

90X 

H 

15(.05) 

14(.05) 

33( . 1 1 ) 

19(.06) 

C 

27S(.93) 

277 (.92) 

266(.89) 

266{.89) 

L 

7(.02) 

9(.03) 

15{.05) 

2.0 

SOX 

H 

11(.04) 

24(.0S) 

50(.17) 

27(.09) 

C 

267(.S9) 

250(.83) 

245(.83) 

240(.80) 

L 

22(.07) 

26(.09) 

5(.02) 

33(.ll) 

90X 

H 

11(.04) 

14(.05) 

28(.09) 

10(.03) 

C 

2S0(.93) 

277(.92) 

271 (.90) 

267(.89) 

L 

9{.03) 

9(.03) 

l(.OO) 

23(.08) 

3.0 

SOX 

H 

23{.0S) 

23(.0S) 

40(. 13) 

19(.06) 

C 

271(.90) 

251{.S4) 

252(.84) 

238(.79) 

L 

6{.02) 

26(.09) 

8(.03) 

43(. 14) 

90X 

H 

S(.03) 

14(05) 

22(.07) 

9(.03) 

C 

2S6(.95) 

277(.92) 

276(.92) 

260(.87) 

L 

6(.02) 

9(.03) 

2(.01) 

31(.10) 

4.0 

SOX 

H 

19(.06) 

23(.0S) 

36(.12) 

15(.05) 

C 

Z71{.90) 

251 (.S4) 

253(.84) 

231 (.77) 

L 

10(.03) 

26(09) 

11(.04) 

54(.18) 

90X 

H 

4(01) 

14(.05) 

21 (.07) 

7(.03) 

C 

2S9( .96) 

277(.92) 

277(.92) 

261 (.87) 

L 

7(.02) 

9(03) 

2(01) 

32( . 1 1 ) 

19 


Table  2 


Average  Length  of  Confidence  Intervals  for  in  P{D>t} 
Exponential  Sojourn  Times 
N=20.  a=0.5.  pj=1.0.  P2=10.0 


Time 

Conf Idence 
Level 

!  BIN 

MLE 

JK 

BT 

0.5 

SOX 

.370 

.144 

.208 

.185 

(.005) 

(.002) 

(.005) 

(.003) 

90X 

.460 

.186 

.271 

.245 

(.007) 

(.002) 

(.006) 

(.005) 

1.0 

SOX 

.545 

.288 

.307 

.280 

(.OOS) 

(.004) 

(.006) 

(.006) 

SOX 

.6S0 

.370 

.400 

.382 

(.OOS) 

(.005) 

(.008) 

(.006) 

1.5 

SOX 

.716 

.431 

.467 

.436 

(.010) 

(.006) 

(Oil) 

(.008) 

SOX 

.SS3 

.553 

.608 

.57S 

(.012) 

(.007) 

(.014) 

(.010) 

2.0 

SOX 

.S15 

.574 

.645 

.587 

(.016) 

(.008) 

(.015) 

(Oil) 

SOX 

1.14 

.737 

.840 

.783 

(.020) 

(.010) 

(.OIS) 

(.015) 

3.0 

SOX 

1.3S 

.860 

1.02 

.832 

(.033) 

(Oil) 

(.024) 

(.018) 

SOX 

1.71 

1.10 

1.33 

1.24 

(.03S) 

(0.15) 

( .032) 

(.025) 

4.0 

SOX 

2.11 

1.15 

1.40 

1.27 

(.06S) 

(.015) 

(.033) 

(.026) 

sox 

2.5S 

1.47 

1.82 

1.70 

(.076) 

(.020) 

(.044) 

(.034) 

Time  Confidence  Coverage  BIN 

Level 


Table  3 


Coverage  Results  for  Confidence  Intervals  for  ^n  P{D>t} 
Exponential  Sojourn  Times 
N=50.  0=0.5,  Pj=1.0.  P2=10.0 


H 

28{ 

C 

245( 

L 

27( 

H 

8( 

C 

282( 

L 

10( 

27 

(.09) 

119(.40)  1 

51 

245 

(.83) 

167(.56)  1 

49 

28 

(.09) 

14(.05) 

0 

27 

(.09) 

70(.23) 

80 

245 

(.82) 

216(.72) 

213 

28 

(.09) 

14( .05) 

7 

9(.03) 

277(.92) 

H(.05) 


H(.05)  22(.07) 

274(.91)  270(.90) 

12(.04)  8^03) 


15(.05) 

268(.89) 

17(.06) 


6(02) 

286(.95) 

8(.03) 


H(.05)  19(.06) 

275(.92)  274(.91) 

11{.04)  7(02) 


10(.03) 

262(.87) 

28(.09) 


‘  ■  • 


Table  4 


Average  Length  of  Confidence  Intervals  for  in  P{D>t} 
Exponential  Sojourn  Times 


Time 

N=50.  0=0 

Conf Idence  BIN 
Level 

5.  Pj=1.0. 

MLE. 

P2=10.0 

JK 

BT 

0.5 

80% 

.216 

.089 

.155 

.127 

(.002) 

(.001) 

(.004) 

(.002) 

90% 

.271 

.114 

.205 

.162 

(.002) 

(.001) 

(.006) 

(.003) 

1.0 

80% 

.319 

.177 

.202 

.181 

(.003) 

(.001) 

(.004) 

(.002) 

90% 

.401 

.227 

.268 

.234 

(.003) 

(.002) 

(.006) 

(.002) 

1.5 

80% 

.416 

.245 

.282 

.262 

(.003) 

(.002) 

(.005) 

(.003) 

90% 

.524 

.340 

.374 

.340 

(.004) 

(.003) 

(.007) 

(.004) 

2.0 

80% 

.524 

.353 

.377 

.352 

(.005) 

(.003) 

(.007) 

(.004) 

90% 

.628 

.453 

.500 

.458 

(.006) 

(.004) 

(.009) 

(.005) 

3.0 

80% 

.750 

.529 

.585 

.545 

(.008) 

(.004) 

(.011) 

(.006) 

90% 

.940 

.579 

.775 

.706 

(.009) 

(.005) 

(.014) 

(.008) 

4.0 

80% 

1.05 

.705 

.799 

.738 

(.015) 

(.006) 

(015) 

(.009) 

90% 

1.31 

.905 

1.06 

.959 

(.019) 

(.007) 

(.020) 

(Oil) 

22 


Table  5 


Coverage  Results  for  Confidence  Intervals  for  Pn  P{D>t}  ^ 

Nonexponential  Sojourn  Times  1 

N=20.  0=0.5.  p,=2.0.  p_=10.0.  p-=2.0  J 

1  ^  J  I 


Time 

Conf idence 
Level 

Coverage 

BIN 

MLE 

(F.xp) 

JK 

BT 

• 

• 

a 

a 

0.5 

80% 

H 

12(.04) 

0(.0) 

105( .35) 

135( .45) 

1 

1 

C 

278(.93) 

53{ . 18) 

183( .61) 

162{ .54) 

i 

1 

L 

10(.03) 

247^82) 

12{ .04) 

3{ .00) 

1 

90% 

H 

O(.OO) 

0{.00) 

73( .24) 

82( .37) 

C 

290(.97) 

109( .36) 

220 { .73) 

217( .72) 

* 

L 

10(.03) 

191 { .64) 

7(.02) 

1(0) 

• 

a 

1.0% 

80% 

H 

16(.05) 

1(0) 

45{ . 15) 

32(. 11) 

1 

C 

254(.85) 

221( .74) 

240( .80) 

241{ .80) 

/ 

L 

30{.10) 

78( .26) 

15(.05) 

27(09) 

j 

i 

90% 

H 

5{.02) 

0(.0) 

29( . 10) 

11(.04) 

j 

■ 

C 

285{.95) 

266(.89) 

263( .88) 

271{.90) 

- 

L 

10{.03) 

34( .  1 1 ) 

8{ .02) 

18^06) 

1 

1.5 

80% 

H 

21{.07) 

o 

00 

40( . 13) 

23(.08) 

C 

260(.87) 

257(.86) 

244( .81) 

232( .77) 

L 

19(.06) 

35(. 12) 

16(.05) 

45(. 15) 

'* 

90% 

H 

5(.02) 

2(.01) 

24( .08) 

7( .02) 

j 

C 

276(.92) 

283( .94) 

273( .91) 

264( .88) 

L 

19(.06) 

I5(.05) 

3(.01) 

29(.10) 

f 

2.0 

80% 

H 

11(.04) 

17(.06) 

39(. 13) 

18(.06) 

t 

C 

270{.90) 

260( .87) 

246(.82) 

232(.77) 

L 

I9{.06) 

23(.08) 

15^05) 

50(.17) 

! 

90% 

H 

11{.04) 

7(.02) 

19(.06) 

7(02) 

C 

284(.95) 

284(.95) 

277(92) 

263{.88) 

j 

L 

5{.02) 

9^03) 

4(01) 

30(. 10) 

3.0 

80% 

H 

13(.04) 

31(10) 

37(. 12) 

I7(.06) 

i 

C 

273(.91) 

254( .85) 

251( .84) 

230( .77) 

r 

/ 

L 

14(.05) 

15(.05) 

12^04) 

53( . 18) 

/ 

J 

90% 

H 

5(.02) 

17(.06) 

19(.06) 

6(02) 

a 

y 

C 

281(.94) 

279 ( .93) 

277( .92) 

261( .87) 

S 

L 

14(.05) 

4(.01) 

4(01) 

33^  1 1 ) 

1 

4.0 

80% 

H 

4(.01) 

44( . 15) 

38( . 13) 

17( .06) 

C 

279(.93) 

245( .82) 

251( .84) 

229^76) 

1 

1 

L 

17(.06) 

11( .03) 

11 (.04) 

54( . 18) 

1 

90% 

H 

4(01) 

23( .08) 

19( .06) 

6( .02) 

3 

C 

296{ .98) 

273( .91) 

277 ( .92) 

260( .87) 

5 

L 

0{.00) 

4(.01) 

4(.01) 

34(11) 

N 

23 


Table  6 


Average  Length  of  Confidence  Intervals  for  in  P{D>t} 
Nonexponential  Sojourn  Times 
N=20.  0=0.5,  Pj=2.0.  P2=10.0  P^=2.Q 


Time 

Conf idence  BIN 

MLE 

JK 

BT 

Level 

(Exp) 

0.5 

SOX 

.289 

.144 

.141 

.  130 

(.003) 

90% 

.359 

.185 

.183 

.175 

(.004) 

1.0 

80% 

.485 

.286 

.269 

.253 

(.007) 

(.093) 

(.006) 

(.004) 

90% 

.604 

.367 

.351 

.331 

(.008) 

(.004) 

(.008) 

(.006) 

1.5 

80% 

.681 

.428 

.451 

.419 

(.010) 

(.005) 

(.010) 

(.007) 

90% 

.850 

.550 

.588 

.552 

(.013) 

(.007) 

(.014) 

(.010) 

2.0 

80% 

.884 

.571 

.643 

.594 

(.013) 

(.007) 

(.015) 

(.017) 

90% 

1.10 

.733 

.838 

.786 

(.017) 

(.009) 

(.019) 

(.015) 

3.0 

80% 

1.46 

.855 

1.03 

.953 

(.034) 

(.010) 

(.024) 

(.018) 

90% 

1.81 

1.10 

1.35 

1.26 

(.043) 

(.013) 

(.032) 

(.025) 

4.0 

80% 

2.34 

1.14 

1.43 

1.31 

(.079) 

(.013) 

(.034) 

(.025) 

90% 

2.84 

1.46 

1.86 

1.74 

(  084) 

(.017) 

(.044) 

(.035) 

24 


Table  7 


Time 

0.5 


1.0 

TTs 

2?0 


Coverage  Results  for  Confidence  Intervals  for  €n  P{D>t} 
Nonexponential  Sojourn  Times 


N=50.  0=0.5. 


P^=2.0. 


P2=10.0. 


P3=2  0 


90X  H  16(.05)  7(02)  20{.07)  12{.04) 

C  274{.91)  282{.94)  269(.90)  264(.88) 

I,  10{.03)  11(04)  11(04)  24(.08) 


3.0  80%  H  25(.08)  34(.ll)  31(.10)  21(.07) 

C  247^  82  )  253(  84  )  248(  83  )  241  (.80) 

L  28(.09)  13(.04)  21(.07)  38( . 13) 
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